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Abstract 

The algebraic approach is employed to formulate N = 2 supersymmetry transforma- 
tions in the context of integrable systems based on loop superalgebras sl(p+l,p), p>l 
with homogeneous gradation. We work with extended integrable hierarchies, which 
contain supersymmetric AKNS and Lund-Regge sectors. 

We derive the one-soliton solution for p = 1 which solves positive and negative 
evolution equations of the N = 2 supersymmetric model. 

1 Introduction 

The well-known examples of the mKdV and sine-Gordon models PP as well as the AKNS 
and the Lund-Regge (Complex sine-Gordon) models |2; can be understood as positive and 
negative flows belonging to the same integrable hierarchy conveniently classified in terms 
of underlying algebraic structure. The basic objects in an algebraic approach to integrable 
models are the loop algebra Q endowed with a grading operator Q and a constant semisimple 
generator E. The grading decomposes the loop algebra into graded subspaces and specifies 
the space of physical fields. 

There already exists an extensive literature devoted to construction of integrable models 
with extended supersymmetry based on superspace formalism, (see e.g. [3 IH El El El El El 
Hm HH II21)- Here, purely algebraic approach without superspace techniques is employed to 
construct a class of integrable hierarchies based on loop superalgebras sl(p + l,p), p > 1 
with homogeneous gradation. Within this framework the symmetry structure of the under- 
lying model is realized as a graded subalgebra of sl(p -|- l,p) defined as a centrahzer /C of 
a semisimple element E {K, = Ker(£^)). The positive and negative grade elements in the 
center of K, generate positive and negative evolution fiows of the extended hierarchy. The 
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positive flows are of tfie generalized AKNS type, while the first negative flow corresponds 
to the Lund-Regge (relativistic) model. By comparing expressions for subsequent positive 
flows we derive the closed expression for the recursion operator TZ. 

The elements in K, not residing in the center give rise to the non-abelian symmetry 
structure of the model. In particular, N = 2 supersymmetry transformations arise in this 
context being generated by specific fermionic generators f± G /C. The supersymmetry algebra 
follows from the fact that the anti-commutator of f± closes on jSj E] 

Furthermore, the algebraic approach allows a systematic construction of Hamiltonians. 
We derive the Lagrangian density for the second flow and the canonical (Dirac) bracket which 
defines the first Poisson bracket structure of the model. Higher brackets can be obtained by 
successive applications of recursion operator. 

We employ the dressing and vertex formalisms to construct soliton solutions. Explicit 
formulas for the one-soliton solution are found for p = 1. The formula is verified to be a 
solution of both the AKNS and Lund-Regge type of evolution equations. 

This work completes the discussion started in [1^], where integrable models based on 
loop superalgebras with principal gradation (like supersymmetric KdV model) were treated 
by algebraic methods. 

Section 2, gives details of algebraic foundation of our construction. We write equations 
of motion governed by the second flow for the s\{p + l,p) model for arbitrary p generalizing 
the results of jHj. Symmetries, conservation laws, recursion operator and the first Poisson 
bracket are presented within the algebraic framework in section 3. In section 4, we deal with 
the first negative flow of the extended hierarchy, which is shown to be invariant under N = 2 
supersymmetry. Next we write explicit the Lund-Regge equations of motion (for the p = 1) 
model using the Gauss decomposition of the zero grade subgroup element B. In section 5, we 
explicit construct one soliton solution by using the dressing and vertex operator techniques. 

2 General Construction of the Model 

Consider the graded Lie algebra s\{p + l,s + 1) composed of bosonic (even) ai, Pt, i = 
1, . . . ,p, t = 1, . . . , s and fermionic (odd), 7, simple roots: 

ai = Ci-Ci+i, i = l,...,p 

A = ft-ft+i, t = l,...,s (1) 

7 = ep+i - /i 

in terms of unit length vectors e^, fa such that ef. = 1 and = —1 with r = 1, . . . ,p + 1 
and a = l,...,s + l. The remaining bosonic roots can be written as 

±(«, + ----f«,_i) =±(e,-eg), ±{Pa + --- + Pb-i) = ±ifa- fb), (2) 

where r,q = 1, ■ ■ ■ ,p + 1 and a, 6 = 1, ■ ■ ■ , s + 1. The remaining fermionic roots have the 
form 

±{ar + --- + ap+i + 7 + A + ■ ■ ■ + A-i) = ±{er - fa). (3) 
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Consider the homogeneous gradation with the grading operator Q = X-^. The affine super 

Lie algebra Q = sl(p + l,s + l) decomposes according to ^ = (BGk with the graded subspaces 
that satisfy [(5,^^] = kQ/. with k E Ij and consist of elements X^^^ = A^X for X G 
sl(p+ + 1). 

The key role in an algebraic approach to integrable models is played by the constant 
grade one generator 

E = E(i) (4) 



obtained by setting n = 1 in the following expression for the semisimple generators of grade 
n 

(P s+l \ 

i=l a=l J 

The centralizer /C = Ker(£') is a subalgebra of Q spanned by the Cartan subalgebra of 
sl(p + 1, s + 1) and the step operators associated to roots 

±{ai + --- + aj) = ±(e, - e^-), ±(/3a + ■ ■ ■ + Pb-i) = ±{fa - ft), 
±{ai H Op + 7 + /5i H h Pa-i) = ±(ej - fa), 

«,J = 1,---,P, a,6 = I,--- ,s + l (5) 

The image I = Im(ii^) is then obtained as a linear combination of step operators of bosonic 
roots 

±(ai H Op) = ±(ei - Cp+i), i = l,---,p 

and fermionic roots 

±(7 + A + ■ ■ ■ + Pa-i) = ±(ep+i -fa), a = 1, ■ ■ ■ , s + 1 , 

respectively. 

The fact that the element E is semisimple insures that ^ = /C © X. In addition, it also 
holds that [J, J] C JC. 

Because supersymmetry requires equal number of bosons and fermions, we impose from 
now on the condition p = s + 1 and denote Op+i = 7, ap+2 = Pi, ctp+s = ^2, - ■ ■ , a2p = Pp-i- 

The super algebra s\{p + l,p) gives rise to a supersymmetric integrable model specified 
by element E of eqn. (@]) and the following Lax operator : 

L = d^ + E + Ao = d, + A., (6) 

where Ax = E + Aq and 

i=l a=l 
i=l a=l 

lies in grade zero sector of X. 
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Now, we propose the zero curvature relation [dt^ + At^,dx + Ax] = 0. We search for 
solution of the form 

=dW+D("-i) + ...D(o), (8) 

where D'^^^ G Gk- The zero curvature relation decomposes grade by grade into the following 
chain of equations, 

= 0, 

■■ ■■ (9) 

The top equation in (jH)) implies that D*^"^ G K, and, consequently, we assume thatL)(") = 

Solving eqns. Q from the highest to the zero grade components we find for t„ = ^2, the 
following equations of motion for the coefficients of Aq: 



dt,h + dlh - 2 ^ [hfhj + FjFj) h = 0, 

p 

dt.bi - dlk + 2J2 {bjbj + FjFj) k = 0, 

p 

+ S^F, - 2 5^ (6,6, + F^F^) F, = 



(10) 



p 

- + 2 ^ (6,6, + = 0, 

i=i 

where i = 1, ■ ■ ■ ,p. These result generalizes equations of motion found in ^Hj in the p = 1 
case. 



3 Symmetries and Conservation Laws 

The symmetries of integrable models will be defined below in a systematic manner in the 
setting of a gauge transformation relating the bare Lax {dx + E) to the dressed Lax operator 
L given by definition (jH} 17 \. 

The dressing operator is defined by 

e = e^(-^'+^'-''+-, 9('>') G g^k 

and satisfies the relation: 

E = Q-^{E + Ao)e + e-^dxQ. (11) 
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Decomposing (fTT|) according to the grading we find the following equations, 



(12) 



Decompose 6^ in its components in the centralizer /C and image X according to 9^ = 
^ic'^ + 4"*^- Then the top eqn. ^ yields 

p 

- f-h -L h _ p -L p p(°) ^ n 

1=1 

The second eqn. (fT^ can be rewritten as 

djt^ + hAo,9^-'^]=0, 

2 (14) 

(—2) 

and can be used to determine a local expression for 9'j^'. 

An element Xm G Ker(i?) of m-th grade generates a symmetry transformation via [T7] : 

(5x„e = (ex^e-i)_ e. (15) 

We will derive the corresponding transformation of Aq using that from the top relation in 
dni) it follows that 

Ao = [9^-'\E] = i[Q,E]),, (16) 

where the right hand side contains projection on the zero grade. 

Employing the definition (fT3j) of the symmetry transformations to the above relation we 
find: 

5x^Ao={[{ex^e-')_^e,E])^ = [{ex^e-')_^,E] 
= {[ex^e-\E])^ = {[[e,E]&~\ex^e-'])^ . 

From equation (fTT|) we find that 

[Q,E]e-^ = Ao + d^eQ-^ 
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and, thus, transformation of Aq can be given by a general gauge transformation formula: 

6x„Ao={[Ao + d,QQ-\eX^Q-'])^ , (17) 
where X„ is an element in /C of n-th grade. For m = 0, we find 

6xoAo = [Ao,Xo]. (18) 

For m = 1, eqn. (fT7|) yields, 

Sx.Ao = [Ao, [9^-'\X,]] + [d.J^-'\X,] . (19) 

In particular, for Xi = E we find that SeAq = d^A^ = dAo/dti, in accordance with the fact 
that the center of /C generates isospectral flows. 
Define two odd elements of Ker(£^) : 

'Mj^m +,(1)77(1) 
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containing grade and 1 generators and Grassmannian parameters e*^°^ and e*-^^ These two 
odd elements satisfy the graded commutation relations: 

[/+,/-] = -2e(°)e«E, 

According to relation (fT?)|l the elements f± G Ker(ii^) give rise the symmetry transformations 
6f^. These transformations satisfy the basic N = 2 supersymmetry relations 

[5u,5f_]Ao = -2e(°)e«9,Ao, [5/±,5/JAo = 0. 

For the particular case of p = 1, i.e. sl(2, 1), we find the supersymmetry transformations 
to be 



5fji 
6f,bi 



jj^ = -e(°)Fi ±6^61 f biFi, 

5f^h = ±e^'^ [dj>i -h ihh + F^h) + h j hFA , 

5f^F, = e^%,Te^^^F, J hF,, 

These are indeed symmetries of the equations of motion (fTUI) . For the sl(p + l,p) case, the 
supersymmetry transformation generated by f± in (j20|l are explicitly listed in the appendix 
A. 

Let us now calculate the Hamiltonian densities. According to fZj, the Hamiltonian 
densities are given by 

^ n— 1 

= - tr = ^ 5^ tr ^A^-k)A(^+k-n)^^ ^ (22) 



fc=0 
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where are defined via 

oo 
k=l 

where dxQQ'^ is the quantity, which entered eqn. (jllll . The symbol tr in expression (j22|) 
denotes a trace for the super algebra sl(p + l,p). 

Now, we explicitly work out the first few expressions. Inserting n = 1 in ()22|1 we obtain: 

n, = -tr(E(°)A(-i)) = ^tr(A^) = ^(6A + F,F,) (23) 

i=l 

where Aq is given by ((Zj). In order to evaluate 7^2 we use (II Hj) to obtain 

p 

n2 = tr (^AodJ^-'^^ = J2 {hdA + FAF.^ (24) 

i=l 

For n = 3, we find 

n, = -^tr(2AoA(-2) + A(-i)') (25) 



where 



and, therefore, 

^ p p p p 

3^3 = Yl + J2 dxFad.Fa + J2 hbjhbi + 2 Yl ^i^^Fb 

j=l a=l i,j=^ i,b=l 

P 

+ (27) 

a,b=l 

3.1 Poisson Structure of the second flow 

The equations of motion (jlUj) can be derived from the Lagrangian density : 
^ p ^ p ^ p ^ p p 

= 2^dtbA--YkdA--YdtFaFa + -Y^-dtFa-Y9^h^^^ 

i=l i=l a=l a=l j=l 

P P P P 

- J2 9xFad,Fa - J2 - 2 5^ lihihFi, - ^ F^FahF^ (28) 

a=l i,j=^ i,b=l a,b=l 

The canonical momenta are given by 

7 



which satisfy the usual canonical Poisson brackets 

{Pb.{x),bj{y)}p = {PiXx),bj{y)}p = -5ij5{x-y), 
{PF^ix),F,iy)}p = {PpSx),F,iy)}p = - y) . 

The momenta dependence on the fields in (j29|) define the primary constraints 

^b, = Pt, - \h = 0, = A, + \h = 0, 
^F.=PF. + \Fa = Q^ $^^ = P^^ + iF„ = 0, (30) 
and accordingly the Poisson bracket algebra is given in terms of the Dirac brackets 

{hi{x),hj{y)}D = -5ij5{x -y), {Fa{x),Fb{y)}D = -5ab5{x ~ y) , (31) 

which together with the total Hamiltonian H3 = J dxTis obtained from (j27jl reproduce the 
equations of motion (fTIH) via formula dX/dt2 = {X, ifsj^- 

A second bracket characterized by P2 can be obtained from the recursive relation 

P2 = 7^Pl , (32) 

where Pi is the first bracket defined in (|HT|) and TZ denotes the recursion operator relating 
consecutive time evolutions, i.e. 



dt^Ao = [E,d,dt^_,Ao] - [Ao,d-'[Ao, [E,dt^_,Ao]]] = nd,_,A 







In components this recursion relation for fiows reads as : 

4 p 

5*™az = EE^S^*'"-^'^'i' ^ = 4, l = l,...,p (33) 
k=i j=i 

where the symbol represents all the fields contained in Aq and introduced in equation 
(fT^ through ^1 ,,■ = bj, C,2,j = bj, = Fj and = Fj. For the case of sl(2, 1) with p = 1 
one finds the pseudo-differential operator expressions for the recursion operator T^j ^ = Tl]'l: 

7^n = d,- 2bd-^b + Fd-^F, TZ^ = -2bd~% 

7^l3 = bd;^F + Fd-^b, 7^l4 = -bd-^F, 

7^2l = 2bd-\ 7^22 = -d, + 2bd-^b + Fd-'F, 

7^23 = bd-^F, 7^24 = Fd-'b + bd-^F, 

7^3l = Fd-^b + bd^^F, 7^32 = Fd-% 

7^33 = -d. + bd-% 7^34 = 0, 

7^4l = -Fd^%, 7^42 = -Fd-^b - bd-^F, 

7^43 = 0, 7^44 = - bd;\ 

where and the pseudo-differential symbol d^^ are acting on all the fields appearing to 
their right. 

In the appendix B we give the general expression for TZ in case of arbitrary sl(p + l,p). 
These expressions allow an easy derivation of the second bracket. 
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4 Relativist ic Model 

Consider the negative grade time evolution equations given by 

dt.,Ao - d,iD^-'^ + + . . . + D^~^^) - [Ao + E, D^~'^ + D^'^^ + ■■■ + D^-^'^] = (34) 
Here, we only consider j = 1. We introduce notation z = t_i and z = —x. Thus, 

dz 

and the above evolution equation becomes 

d,Ao + BD^-'^ - [Ao + E, = , 

with the following solution 

^(-1) = Ao = dBB-^ (35) 

where B is the group element of sl(p+ l,p) satisfying the Leznov-Saveliev's equations 

d{B-^d,B) + [E'--^\ B-^EB] = 0, d,{dB B~^) + [BE^-^^B~\ = . (36) 

It is customary to think about coordinates z as light-cone coordinates. Such association 
allows us to interpret the model as being relativistic. 

The group element B transforms under the symmetry transformations as 

6x^B = - (eX„0-i)^ 5 , n > . (37) 

For n = and n = 1 the above relation takes form 

6x,B = -XoB, 6x,B = [Xi,e^-'^]B . (38) 

In particular, for Xi = E the above definition yields 

SeB = - {QEQ-^)^ B = -[e^-^\E]B = -A^B 

recalling that Aq = —d^B B~^ we obtain 5eB = d^B, confirming again that E generates the 
first positive isospectral flow ti = x. 

These transformations leave the above equations of motion ()36p invariant. In particular, 
the equations of motion ()36|) remain unchanged under the N = 2 supersymmetry transfor- 
mations Sf^. 

Notice that eqns. yields chiral currents associated to the Kernel of E, i.e. /C = {X e 
sl(j9 -|- l,p), [X, E] = 0}. In order to make connection with the super AKNS model ((Tj), i.e. 
Ao € X in (jH^j) . we impose the following subsidiary constraints 

tr {XdBB-^) = tr {XB-^d,B) =0, X e IC (39) 
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4.1 Relativistic Super Lund-Regge Model 

In order to define the relativistic member of tfie so-called super Lund-Regge hierarchy we 
consider the zero grade group element B E Gq = SL{2, 1) parametrized as 

where we are using the basis E±ai, E±a2, E±(_ai+a2)y " H and 0^2 ■ H. Changing to the 
natural variables 

~ 'Pl+'P2 _ _V2 ~ „ 'P2 

tjj = ^e 2 , gi= gie 2 , f^ = f^e ^ 

X = Xe ~, g2=g2e~~, f2 = f2e~~ (41) 
we find for the currents 

H (ai + 02) ■ H 

Ja2 H ■ ~l~ JaiE—cti ~l~ Ja2^—a2 ~^ ^ai+a2 (cn+a2) ' 
J = ""^ = J-aiEai + 1-02^02 + J-ai-a2Eai+a2 + ■J{ai+a2)-H («! + ^2) " -f^ 

— Ja2 H Ot2 ' H -\- Ja^E^ai + Ja2-^-a2 + <^ai+a2-^-(ai+a2) (^2) 

with the coefficients J-^i, . . ., Jai+a2 shown in Appendix D. 
The constraints 

J±{ai+a2) — J±{a-i+a2) — J{ai+a2)-H = J{ai+a2)-H = Ja2H = Ja2H = (43) 

lead to the following equations for the non local fields 

dzfi = ^fidz'^2 + g2[dzX - ^xidz^pi + dz(p2)], 

dz9i = i^dj2 + ^9idzip2-^4^f2dz<^i, 

1 - I - 
dfi = xdg2 + -fid<f2 - -X92d(fi, 

dgi = ^gid<^2 + f2[d^ - ^'ipid^i + 9:^2)], 



dz^2 

dipi 



^[dzX{l + 92f2) + \X92dzf2] 

i + ^x(i + 1^72/2) 

i^dzXi'^ + 1^2/2) - g2dzf2 - lipX92dzf2 

i + ^x(i + 1^72/2) 
X[a^(l + ^2/2) + 1^^2/2] 

l+^X(l + k2/2) ' 



. ^ X^^(l + 1^72/2) + (i^X+l)/2a^2 

i + V^x(i + W2) ^ ^ 
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The equations of motion are given in the Leznov-Savehev form ()36p. after ehminating the 
non-local fields with help of equation ()44|). In components they read 



1- 1- - 1-- 

(1 + 92f2)[dd^x - ^dx{d,^i + 8,^2) - ^x{dd,^i + dd^^2) + -^{d^i + 9^2) 

{d.X - \x{d.Vi + dz^2))] + d{g2f2)[d.x - \x{d.Vi + ^.^2)] + X = 0, (45) 

ddzg2 + ]^dg2dz(pi + ]^g2ddz(pi - ^dzg2dipi - ^g2dipidzipi + (1 + 'ipx)92 = (46) 

(1 + g2f2)[dzdip - ^d^tpidipi + dip2) - ^ij{dzd(fi + 8^8(^2) + ^(^^V^i + d,ip2) 

{d^ - l^dVi + dy^2))] + dM2)[d^P - ^HdVi + 0^2)] + = (47) 

dzdf2 + ]^dj2d^i + ]^f2dzd^i - ]^df2dz^i - ^f2dz^id^i + (1 + V^x)/2 = (48) 



where the non local fields and (p2 are given in equation 

4.2 Connection between AKNS and Lund-Regge Variables 

On basis of (jHSj) we find the following non-local relations between the AKNS and the rela- 
tivistic Lund-Regge variables : 

hi = J_„i = — — — dip - -ipidifii + difi2) 
Fi = J-0.2 = e-^'^^ (df2 + \f2dip1 



bi = Ja^ = e-^^^'^^+^''^ (^dx+lx{dy^i + d^2)-xf2dg2-^xdvi92f2 + e-^^^ 

Fi = Ja2 = ei^^ (^a(72-^^?29<^i + e-^('^^+'^=^)/iJ_«,^ (49) 

Notice that they are all non local fields in terms of the Lund-Regge variables due to the 
presence of <fi,<f2 and /i defined by (jUj). Also, we have seen that in terms of the AKNS 
variables, the models are invariant under supersymmetry transformations ()21|). 

5 Soliton Solutions 

5.1 Dressing Transformations and Vertex Operators 

The dressing transformation and the vertex operators method represents a powerful tool 
for the construction of solitons solutions of integrable models. The dressing transformation 
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relates two solutions of the equations of motion written in the zero curvature representation. 
In particular, it relates the vacuum and the 1-soliton solutions by a gauge transformation, 

A, = Q^Al^'Ql^ - (9^e±) Ql^ (50) 

where 

0_ = e*(-i)+-, 6+ = e^(°)e^(i) • • • (51) 

The zero curvature representation implies for pure gauge solutions: 

j^vac ^ ^o,ToT,-\ = -d,TT~' (52) 

which allows the following relation 

T = e±To, I.e. e+T = e_Tof7, (53) 

where G G is an arbitrary constant element of the corresponding affine group. Suppose Tq 
represents the vacuum solution. 

To = exp(-t„E(")) exp(-xE), n > 1 (54) 

i.e.. 

As consequence of ()50|) with ()55|) and (jM]) we can determine G±. Consider for instance eqn. 
()50|) for Arc and G_. It determines its zero grade component t{—l) through 

Ao = [t(-l),E] (56) 

Notice, that if we enlarge the loop algebra by adding central terms, consistency requires 
the introduction of a additional field u associated with the extension Aq ^ Aq + dxvc in 
constructing the Lax operator. For the specific sl(2, 1) model: 

Ao = hE^^; + WE^X + Fi^W + F,eX + d,vc (57) 

we find 

= + i^^E^) + _ F.E^l + d,v]^{a, + a^) ■ if^'^) (58) 

From we find the solution 

To = =< Xo\TogT^^\Xo >, 
n = b2To =< XoIe'lI^ToqTq^IXo >, 
T2 = b2ro = -<Xo\E^^^TogT,-'\Xo>, 
rg = F,To=<Xo\E%TogT^'\Xo>, 
n = Firo =< Aol^g^ro^^To-^IAo >, 

(59) 
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For the relativistic model let us consider eqn. ()50p for and 6^. We find 

Ao = -9,e"(o) e-"(°) (60) 
Since,for the affine case, Aq = —d^BB^^ — d^vc, 

From ()53p we find the solution 

< K\Be^'\\t >=< Xr\TogT,-'\Xt > (62) 
We now chose |Ar >, |Af > to be the following weight states 

1 1 
|p>=|-ai>, \e>=\- -ai-a2> (63) 

and obtain 

<p\Be'''\p> = roe5('^^+^^) =< p|To(7To-V >, 
< elBe^^le > = roe^^'^i-'^^) e\TogTo^\e >, 

<p\Be''Eo.,\p> = -^~g2Toe"^^^^+^''> =< p\TogT,~'E^,\p>, 

< p\Be-'E_^, \p> = roe5('^i+^2) + 1^2^i^ =< p\TogT^^E^^, \p >, 

< plE^^.Ee'^^lp > = -i/2roe^('^^+'^^) =< p|E_,,To^7To-V >, 

< p\E^,Be'''\p > = roe^('^i+'^2) (^i, + i^2^i^ =< plK^To^^To-^ >, 



(64) 



5.2 Vertex Operators and Solutions 

In order to evaluate the matrix elements in the r.h.s. of eqns. ()59|) and ()64|) it is instructive 
to introduce the eigenvectors (^^(7)) of E"^, i.e., 

[E^"'\Hl)] = frh)H7), (65) 

such that we can classify the soliton solutions in terms of the constant group element g = 
g-^i(7i)g-7^2{72) . . . g^j^^ hence 

TogTo^ = exp (pi (71)^1 (71) + P2(72)^2(72) + ■ ■ ■) (66) 
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where ^^(7^) = exp (-tn/f - xfl). 

For Q = sl(2, 1) model, E^^^ = (ai + 0:2) ■ H^^'^ and its eigenstates are given by 



— Q2 



J-i(7) = 5^(i?£)+aiE 

-^2(7) = E(^-t+«2i?£ 



7" 
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(67) 



with eigenvalues /|"(7) = 7™, fTil) = ~7 We now consider soliton solutions associ- 
ated to g = exp(jFi(7i)) exp (^2 (72)) and 

TogT,' = (1 + pi(7i)-Fi(7i)) (1 + P2(72)^2(72)) (68) 

leading to the following matrix elements 

To = <0\TogT,-'\0>=l + TM, 

n = <0\Ei\TogT,-'\0>=-^,p',, 

T2 = <0|To^?To-iE(;i)|0>=7ip'i, 

ra = <0\E^}lTogTo'\0>=-a2i2P2, 

n = <0\TogT,-'Ei-'^\0>=aijip[ (69) 

where 6 = 0102, Fq = (^^l^^^p and F = (1 — 6)Fo. We therefore find the solution for the non 
relativistic model (jHH), 

7iP'i , 72P2 



h 

El 



i + rplp'^' 
72P2 



bi 



-02 



1 + Fp>'2 ' 
7iP'i 



l + Fplp'^' 



El = ai 



l + Fpip^ 

P'2(7l) = 



(70) 



3*272+722; 



where for the non relativistic model p'i(7i) = e 

For the relativistic model described by equations of motion pH|l - ()48|l we find the following 
matrix elements 



< p\TogTo V > = 1 + P1P2 

< elTogT^'le > = 1 + pip2 



7i 



(71 - 72)^ 
7i 

(71 - 72)^ 



71 - ^(71 + 72)^ 
71 + ^(71 - 372)^ 



< p\TogT,'E^'l 

< p\E^STogTo' 

< p\TogT,-'Ej^l 

< plE^XngTo' 

< p\TogT,"'E%_^^ 

< p\E!^;^^jogTo' 



P> = Pi, 

P> = P2, 
1 

2' 
1 



P> 
P> 
P> 



-^(^IPl, 



(71) 



-2«2P2, 

1 7i 

o«lPlP2 



(71 - 72) ' 

7i 
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which allow us to determine the solution 



1 + jro (i 



^ (71+72) 

To 



P1P2 



_ ^ + ^° + ^^^) ^^^^ _ , _ b (72 - 7l)roPlP2 

&71P1P2 



^0 2 72(1 + roPiP2) 



X 




2(7i-72)(l + 5ropiP2);' (72) 



&71P1P2 



2(7i-72)(l + 5roPiP2) 



71P1P2 , _ ^ 71P1P2 



fi'i = a27 h = ai 

(71-72)7-0 (71-72)^0 

^2 = ai^e-i^S /2 = a2^e-^^^ 

From the solutions given in ()72p it is easy to verify the consistency of the constraints (jH^ 
(or (HH)). Let us take for instance 

dfi = xdg2 + ^fidf2 - ^g2xd(pi (73) 

Substituting (f?^ in both sides of eqn. (fTSj) and identifying t_i = z, x = —z using the 
second eqn. (f?^ together with the fact that aib = a2& = we find that the constraint (f73j) 
is satisfied. The same can be shown for all eqns. ( i.e. all constraints are satisfied). 

Finally we verify the connection between the non-relativistic and the relativistic variables 
given by eqns. (jl^ . Substituting the solutions (f?^ in (jl^ and comparing with (f7n|) we 
conclude that they agree if we associate 



7% - = ^, x^--z (74) 
7 7 



which implies p(7) — > p'{'~f). 
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A Supersymmetry transformation for sl(;) + 1,^) 

Let f± be the supersymmetry generators from equation (j20|) . Using the symmetry transfor- 
mations given in equations (|T^ and (fT!?| and 

i=l 



and 



a=l 



p p 



j=i, i=i j=i 
p p p p 

j=l a=l j=l a=l 

p p p 

a=l, aj^b b=l a=l 

obtained from equation (fT^ we find 

Sf^Fr = ± (^d,br - ^E / (^^-^^^ + ^^■^-) j + / j 

= hre^^^ T^Fj h^F^e^^^ (77) 

B Recursion Operator for sl(^ + 1,^) 

Using the same notation as in equation (jHHj) we write down expression for the recursion 
pseudo-differential operator Ti!'^^. vahd for a general case of sl(p+l,p). Here, the /, j = 1, . . .,p 
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and indices i,k — 1, 2, 3, 4 label the four modes bj, bj, Fj, Fj. 

= 5i, [d, - J2 hd-% + Fbd-'F^j - bid-%, 

n'^2 = -hVbi-bid-'b„ 
ng = hd-'F, + F,d-\, 
nl = -kd;'F„ 
ng^b,d-% + bid-%, 

= 6i, (^-d, + bkd-% + Y Fbd-'F^ + bid-'b„ 

ng = bid;'Fj, 

T^'L = FA\ + bid-'F^, 

nl = Fid-\ + b^d-'Fi, 

n'i, = Fid-'F, - F,d-'Fi, 
T^l = -d-'hFu 

n% = -Fid-'b,-b,d-'Fi.), 

n'l, = Sij (d, -J2hd;% + ) - Fid^'Fj. 

\ k b J 

C Case of S(3, 2) 

Consider the sl(3, 2) case where 

+ + F.E^Z + F,E^:i^^ + F.E^X-o.. 

and propose the zero curvature representation 
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from where we 



find the following solutions for equations 



= (D(2))^^^ = (ei + e2-/i-/2)-i/('\ 

+F,ES + F.E^^l,^ + F,E% + F^E^X_^^ 
{D^'^)iM = -d,hE^\_^^ - d,hE% + dME^^l^^ + dJj^E^S 

- 9..F2<V., + d^F,E%^ + d^F,E^X_^^ 
(DW)xe. = -{hihi)hf^ - ihh + hh)hf^ + {F,F^ + F^F^)hf + {F^h)hf - {h,h)Et 

-{biF,)E 

—ai— 02—03 

(62&i)<) - {b2F,)E 

—02— 03 

-(&2-p2)-Ei°a2-03-04 ^ (^l-^l)-^i°i+02+03 - (&l-^2)-Ei°Vo2+03+04 ~" (^2-Fl)-Ei1+c 

-(^2F2)4t„3+„, - {F^F2)E^X - {F,F,)E^S 



dtM + 0% - 2 (6161 + 62^2 + i^iA + F2F2) &i = 
&2 + <9^&2 - 2 (6161 + 62^2 + F^Fi + F2F2) &2 = 
dtj}i - 9^61 + 2 (6161 + 62^2 + FiFi + F2F2) &i = 
5*2^2 - d% + 2 (6161 + 62^2 + FiFi + F2F2) &2 = 

^i^Fi + ^^Fi - 2 (6161 + 62^2 + i^2i^2) i^i = 

^^2^2 + 92^2 - 2 (6161 + 62^2 + F^Fi) F2 = 

9^2^! - dlF^ + 2 (6161 + 62^2 + i^2i^2) i^i = 

dt.h - dlF2 + 2 (6161 + 62^2 + FiFi) F2 = (81) 



(80) 



and the equations of motion 



Consider 




(82) 



SxoFi 
SX0F2 
SxoFi 
SX0F2 




(83) 
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Of course when bi = bi = and F2 = F2 = we recover the p = 1 or SL{2, 1) case proposed 
in ref. ^14J corresponding to Aq = 62-E'ai + &2-E'-ai + FiE^^ + FiE_a2, i.e., 

dtM + dlb2 - 2 (62&2 + FiF,) 62 = 
5t2&2 - 5^&2 + 2 (62^2 + FiFi) 62 = 
Fi + dlF, - 26262 Fi = 



dt,F, - dlF, + 2b2b2Fi = 



(84) 



D SL(2,1) super-currents 



The components of the sl(2, 1) currents J = B ^dB and J = dBB ^ given in equation (jl^ 
read in terms of the fields defined in (jlOj) and (PT|) : 

= e^'^^(^9/i-^/i9<^2-^72e-^(^^+^^)j„,^, 



2 



j^^^ ^ ^ _^ ^^^^^^ _ l^0v^g^f2 + giei'^'Ja 
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and 

Ja,.H = dv2 + dg2f2 + lf2g2d^i-xe-'^'^'''^'''^J-a„ 
Ja, = e-5('^i+'^^) ("sx + ^x(a</^i + 8^,2)- Xf2dg2 - ^X^</^i^?2/2 + fie"^^' J-a, 



(85) 



References 



[1] A. Chodos, Phys. Rev. D 212818 (1980); C. A. Tracy and H. Widow, Commun. Math. 
Phys. 179 1 (1996) solv-int/9506006 ; L.A. Ferreira, J. L. Miramontes and J. S. Guillen, 
Journ. of Math. Phys. 38 882 (1997) hep-th/9606066 ; J. Dorfmeister, H. Gradl and 
J. Szmigielski, Acta Applicandae Math. 53 (1998) 1; D. Fioravanti and M. Stanishkov, 
Nucl. Phys. B 591, 685 (2000) |hep-th /0005158 ; 

[2] H. Aratyn, L. A. Ferreira, J. F. Gomes and A. H. Zimerman, J. Phys. A 33, L331 
(2000) nlin.si/0007002|; A.M. Kamchatnov and M.V. Pavlov, Phys.Lett. A 301 269 
(2002) nlin.si/0208025j 

[3] F. Delduc and L. Gallot, Jour, of Math. Phys. 39, 4729 (1998) |solv-int/98020l3|. 

[4] J. O. Madsen and J. L. Miramontes, Commun. Math. Phys. 217, 249 (2001) 
|hep-th/9905103|. 



[5] Z. Popowicz, J. Phys. A 29, 1281 (1996) |hep-th/9510185 



[6] A. S. Sorin and P. H. M. Kersten, Lett. Math. Phys. 60, 135 (2002) |nlin.si/0201026| . 

[7] F. Delduc and A. S. Sorin, Proc. of 2002 Workshop on Integrable theories, solitons and 
duality, eds. L.A. Ferreira, J.F. Gomes and A.H. Zimerman , J. High Energy, PRHEP- 
unesp2002/044, ^nlin.si/0206037| 



20 



[8] E. Nissimov and S. Pacheva, Jour, of Math. Phys. 43 2547 (2002) |nlin.si/0103055l 



[9] E. Ivanov, S. Krivonos and F. Toppan, Mod. Phys. Lett. A 14, 2673 (1999) 
|solv-int/9912003 . 

[10] F. Delduc, L. Gallot and A. Sorin, Nucl. Phys. B 558, 545 (1999) |solv-int/ 9907004|. 

[11] F. Delduc and L. Gallot, Commun. Math. Phys. 190, 395 (1997) |solv-int/ 9609 0081. 

[12] E. Ivanov and S. Krivonos, Phys. Lett. A 231, 75 (1997) hep-th/9609191|. 

[13] H. Aratyn and C. Rasinariu, Phys. Lett. B 391 99 (1997) | hep^h/9 608107|; H. Aratyn, 
A. Das and C. Rasinariu, Mod. Phys. Lett. A 12 2623 (1997) |solv-int/9704lT9| 

[14] H. Aratyn, A. Das, C. Rasinariu and A.H. Zimerman, in Supersymmetry and Integrable 
Models" , H. Aratyn et.al. (Eds.), Springer- Verlag, 1998 (Lecture Notes in Physics 502) 

[15] H. Aratyn, J. F. Gomes and A. H. Zimerman, Nucl. Phys. B 676, 537 (2004) 
|hep-th/0309099|; hep-th/04 09171j ; Contribution to the 11-th International Confer- 
ence on Symmetry Methods in Physics (SYMPHYS-11), Prague, Czech Republic, 
|hep-th /040823l| 

[16] T. Inami and H. Kanno, Commun. Math. Phys. 136 519 (1991) 

[17] H. Aratyn, J.F. Gomes, E. Nissimov, S. Pacheva and A.H. Zimerman, ^'"Symmetry Flows, 
Conservation Laws and Dressing Approach to Integrable Models" in Integrable Hierar- 
chies and Modern Physical Theories, Ed. H. Aratyn and A. Sorin, Kluwer Academic 
Pubhshers (2001), P. 243-275 



21 



